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SUMMARY OF CONTENTS 

PART I - What the DA-1 Is and Does Pages 5 - 10 

Accessory DA-1 permits the direct programming 
of the solution of differential equations on the Bendix 
G-15 Computer, Characteristics of the DA-1 and 
the method of its use are summarized. 

PART II - How the DA-1 Works Pages 11 - 14 

The method of operation utilizes the "Integrator" 
concept, familiar to users of analog computers. 
In this case, the integrators are numerical regis- 
ters on the drum of the G-15 Computer. The tech- 
nique andprinciples involved are described in some 
detail. 

PART III - What Integrators Can Dp Pages 15 - 30 

The integrators in the DA-1 are versatile. In addi- 
tion to the integration operation, they may be made 
to act as decision elements, simulate servo action, 
or add variables together. By their interconnection, 
a great many mathematical functions can be gener- 
ated. The ways in which integrators may be used 
are described with many examples of function gen- 
eration. 

PAR T IV - How to Map Solutions Pages 31 - 43 

A "map" is a flow diagram of interconnections be- 
tween integrators and represents the solution of a 
differential equation, or set of equations. The 
preparation of a map is the first step in program- 
ming. A typical differential equation is mapped 
three different ways in the text in order to illus- 
trate the techniques involved. 

PART V - How to Scale Numerical Values Pages 44 - 48 

"Scaling" refers to the selection or determination 
of maximum values assumed by integrands and of 
numerical values of incremental inputs to integra- 
tors. The scaling of the flow diagram is the second 
step in programming. The methods and rules for 
scaling are specified and explained. 



PART VI - How to Code and Obtain Numerical Solutions Pages 49 - 67 

The third and final step in programming is putting 
the scaled flow diagram in codedform for insertion 
into the computer. Instructions are set forth for 
coding, for obtaining readout of the solution, and 
for operating the computer. 



APPENDIX 



Pages 68 - 71 



No manual is complete without an appendix. Des- 
cribed in this one is the physical form of informa- 
tion storage on the G-15 memory drum when Ac- 
cessory DA-1 is used. Since the modification of 
information may be programmed when its form is 
known, programs which simultaneously take advan- 
tage of the general computation and differential an- 
alyzer facilities of the G-15 system may now be 
written. 
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Plot of an elliptical spiral generated with accessory 
DA-1. Four integrators were required: two to gen- 
erate the curve, and two to control the graph -plotter. 
The simple program is described in the text. 



PART I - WHAT THE DA-1 IS AND DOES 

Purpose of The Digital Differential Analyzer Accessory, Model 

Accessory DA-1, permits the simple programming of the numer- 

ical solution of problems which can be expressed in the 
form of ordinary differential equations. These include 
the solutions of linear and non-linear ordinary differ- 
ential equations, linear and non-linear simultaneous 
equations, solutions for roots of transcendental equa- 
tions, and the simulation of real systems. 

In order to solve differential equations on a general 
purpose computer without a differential analyzer attach- 
ment, it is necessary to transform the equation into 
suitable arithmetic form by use of techniques of numer- 
ical analysis and then to formulate coded instructions 
for the computer. The programming time required is 
extensive and the programming knowledge required is 
specialized and often sophisticated. By use of Acces- 
sory DA-1, the differential equations, in a slightly 
altered form, may be placed in the computer andtheir 
digital solutions realized. Since the solutionis digital, 
much greater accuracy may be obtained than is possi- 
ble from earlier-designed differential analyzers which 
are of an analog nature. 



Background The basic principles of a differential analyzer for the 

solution of differential equations are not new - they 
were conceived by Lord Kelvin 100 years ago, des- 
cribed by Sir William Thomson 80 years ago, and put 
into practical mechanical form by Dr. Vannevar Bush 
30 years ago. But the means of execution embodied in 
the DA-1, the use of digital techniques and the associa- 
tion with a general purpose computer, permit an accu- 
racy, versatility, and ease of use which greatly extend 
the range of usefulness of this type of machine. 

The basic computing unit in the differential analyzer is 
the "Integrator". An integrator may be thought of as a 
physical unit (electronic digital registers in the DA-1; 
mechanical shaft assemblies in the early Bush machine) 
which has two input lines and one outputline. The units 
may be interconnected so that the output lines of some 
integrators become the input lines of other integrators. 
The inputs are differentials of dependent and independent 
variables arising from the differential equation being 
solved; in this manual they are called dy and dx, res- 
pectively. The output of the integrator is a differen- 
tial, dz, so that 

dz = y dx 



The value y, which is the accumulated value of the dy 
inputs, is held in a register within the integrator. 
Therefore, if the output line of the integrator is 
connected to the dependent variable input line of a 
second integrator, we may obtain in the second inte- 
grator 

z = / dz = f ydx 

Since differential equations usually involve constants, 
the output of an integrator may be multiplied by any 
arbitrarily chosen constant value before being used. 
We may, then, obtain with two integrators 

dz = K ydx (Equation 1) 

z = K /ydx (Equation 2) 

By programming interconnection of integrators of 
this nature and by selection of initial values for the 
"y" term in each, any ordinary differential equation 
may be solved. 



Method of 
Solution 



A frequently-used approach is to isolate the highest- 
order derivative in terms of the other variables. 
The existence of this derivative in an integrator is 
assumed and by repeated integration is used to pro- 
vide all the lower-order derivatives. The lower - 
order derivatives are then combined according to 
the conditions of the equation to find the highest-order 
derivative. The resulting closed loop may then be 
put in the computer in coded form and solved. An 
illustration of the technique is given in the example 
below. 



Example 



The equation to be solved is: 

!LZ + . 297 uv 2 p- - u Z v + sin u = 
, 2 du 

du 



Isolating the highest order derivative: 



d 2 v 



du 



2 dv 2 

2Q7 uv -r- + u v - sin u 
du 



By integrating once with respect to u: 
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This equation may be put in differential form; 

d~ = -. 297 d /uv 2 dv + d /u v du + d cos u (Eq. 3) 

Since theinputsand outputs of integrators are differenr 
tials, Equation 3 now describes a relationship between 
integrators. A machine independent variable called dt 
is available in the DA-1. This variable provides du. 

4^ is assumed to exist in the y register of one integra- 

dv 
tor. By development from -3— and du, the right-hand 

terms in the equation may be obtained as the outputs 
of integrators. These outputs may then be combined 

in the manner of Equation 3 to form d-r-. By making 

d^ the dependent differential input into the integrator 

in which -r^ was assumed to exist, the feed-back loop 
du 

is closed and the layout of the solution is completed. 

The procedure ia summarized in block form in Figure 
2. Each block represents one integrator and each 
functions in the manner of Equations 1 and 2. The 
solution of this problem is described in detail on page 
31. 

The diagram of the solution is then put in coded form 
and inserted into the computer. By use of a standard 
program preparation routine, the initial values of the 
■y" terms in the integrators and the values of the con- 
stant multipliers may be expressed as decimal num- 
bers. The interconnection paths of the integrators 
are specified numerically in the program code. 

Since Accessory DA-1 provides 108 integrators and 
108 constant multipliers (a much greater number than 
heretofore available in digital differential analyzers), 
it may be used to obtain numerical solutions of com- 
plex problems involving high-order differential equa- 
tions or sets of simultaneous equations. 



Character- Distinctive features of the DA-1 are listed below. 

istics 

108 integrators and 108 constant multipliers are pro- 
vided. 

The primary functions of an integrator are those des- 
cribed in Equations 1 and 2; but, an integrator may 
also be used as a servo, as a decision element, or for 
the addition of variables. 



The outputs of all integrators are available to be used 
as inputs by every integrator. Interconnection is by 
numerical coding and not by patchboard. 

During a single step of computation, an integrator 
receives dy and dx inputs and accumulates both dy and 
ydx products. Whenever the sum of ydx products ex- 
ceeds a pre-assigned quantity, the integrator has a dz 
output. Accuracy and speed of computation has been 
increased by making the dz output ternary rather than 
binary: dz may assume a value of + one increment, 0, 
or -one increment. 

An integrator may be coded to receive a dx input from 
one other integrator; it may be coded to receive dy in- 
puts from any number of other integrators. The value 
of the dy input during one step of computation must 
range between -7 and +7 increments. 

The set of 108 integrators is processed 34 times per 
second. The speed of computation is independent of 
the number of integrators used in the program. The 
speed of solution of a problem is inversely propor- 
tional to the degree of precision required. 

The general computing facilities of the G-15 maybe 
used in conjunction with the DA- 1 for the programming 
of highly involved problems. 

Empirical functions may be used. 

A standard program preparation routine provides the 
following input-output characteristics: 

Initial values of integrands and values of constant 
multipliers are expressed in decimal form as positive 
or negative floating decimal point numbers. Precision 
of computation remains optional. Up to 7 decimal 
digits of precision is available. Exponent range is 
from 10" 4 to 10+9. Computed results are typed out 
as floating point decimal numbers. 



Physical and The DA-1 is controlled from the G-15 General Pur- 

Electrical pose Computer. The DA-1 has no controls of its own 

Specifications and may not be used independently. 



Input can be programmed from normal G-15 equip- 
ment. This includes typewriter, punched paper tape, 
punched cards and magnetic tape. Empirical function 
input may be from any G-15 input device. Output may 
be tabulated by normal G-15 equipment or plotted 
graphically. 

The graph-plotter plots the relationships between any 
two variables generated by the computer in 0.01 inch 
increments on a standard, sprocketed roll of paper, 
1 foot by 100 feet. It is an optional accessory to the 
DA-1. 

The DA-1 is 22 inches deep by 24 inches wide and is 
60 inches high. It weighs 300 lbs. Power input is 1 
kva, 110-120 volts, 60 cycles, single phase. Cooling 
is by internal forced air. 
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PART H - HOW THE DA-1 WORKS 



It has been stated that numerical integration can be 
performed with two integrators. An incremental 
output, dz, (which is equal to either +1, - 1 or zero) 
is obtained from one integrator so that 

dz = ydx 

and the incremental outputs are accumulated in a 
register in a second integrator in order to obtain 



z - / dz 



ydx. 



The procedure used may be explained graphically: 






Y 



/f\ 



y* 



dx 



a#£X *&-f K,</x 



Figure 3 

Referring to Figure 3, consider that we are evaluat- 
ing the definite integral of Y m with respect to x, by 
continuously accumulating the area under the curve 
Y m = f(x) as x increases. 
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Initially, the curve is at point y Q , x Q . The curve is 
approximated by incremental changes in y and x. The 
incremental changes are called dy and dx on the dia- 
gram. If Y m is the instantaneous value of y, each 
time x increases by one dx increment, a rectangle of 
width dx and of height Y m will be added to the accumu- 
lated area. The incremental changes in Y m »which may 
consist of from one to several increments at a time, 
will be rounded off in a manner which reduces the tri- 
angular error made at the tops of the rectangles. The 
method of this round off will be discussed later. Al- 
though we must keep track of the value of Y m in order 
to know how high a rectangle to add with each dx change, 
we need not keep track of the entire accumulated area. 
We may simply send out notification each time an in- 
crement of area has been covered and allow the total 
area to accumulate elsewhere. It will then be neces- 
sary to keep track of only Y m and of any fractional 
part of an increment of area. In Figure 3 we will con- 
sider that a rectangle of width dx and of height +1 (the 
full height of the diagram) will be our increment of 
area represented by a dz pulse. 

It is in this manner that the DA-1 integrator operates. 
Each integrator, which is represented by a block as 
shown in Figure 4, receives two incoming lines of in- 
formation. One of these, the dx line, notifies the in- 
tegrator each time that x has increased or decreased 
one increment. The dy line notifies the integrator of 




- dz 






Figure 4 

the corresponding changes in Y m . The output of the 
integrator is the dz line through which the integrator 
signals each time an increment of areahas been trav- 
ersed. If the dz line is fed into another integrator as 
its dy input, this second integrator will then contain 
as its y value the total instantaneous value of the area Z. 
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It is necessary for the integrator to contain two regis- 
ters (see Figure 5). One, called the Y register, al- 
gebraically accumulates the value of Y m . The other, 
called the R register, contains the fractional part of 



1 5 m 



i 



dz 



R /?£<5/S r£& 



Figure 5 

increments of area; the value in the R register will 
range between and one increment. 

The actual digital operation is accomplished in a 
straight -forward manner, dy increments are added to 
Y m algebraically. Each time a dx pulse is received, 
Ym is added into the number which exists in the R reg- 
ister. Since neither register is cleared in the process, 
the R register will periodically overflow. It is the 
overflow of the R register which causes the dz output 
signals. Each dz pulse corresponds to the increment 
of area used in the graphical illustration. 

The multiplication of an integrator output by a constant, 
k, is similarly done by two registers. One, the K 
register, has no modifying inputs . The other, the 
K r register, has the value in K added to it whenever 
a dz output occurs. The overflow of the K r register 
provides Kdz output signals. 

Electronically the Y, R, K and K r registers are par- 
allel memory channels on the magnetic drum of the 
G-15 computer. Each word-time in the channel(s) 
represents one integrator. Since there are 108 word- 
times in a drum channel, we may store information 
for 108 integrators on the drum. 

Consequently, integrators are processed serially, be- 
ginning with integrator 00. A single processing of all 
108 integrators, in which dz outputs are transferred 
between integrators, is called an iteration. 
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As a second consequence, the time required to compute 
a solutionis not a function of the number of integrators 
used in the program, but is a function of how soon in- 
cremental inputs to integrators cause incremental 
outputs. The amount of variation in the y value, or 
integrand, to be represented by a single incremental 
input is optional. Naturally, the smaller this incre- 
mented variation in the integrand is made, the greater 
the precision of the solution. Each additional binary 
digit of precision doubles the time required to compute 
a solution. 



Round-off As has been stated, a round-off correction is made in 

Correction each integrator. The correction consists of starting 

with one^half full value in each R register. Now a dz 
overflow will occur from the addition or subtraction of 
more than one -half full value into the R register. 

Therefore, if dz is used as the dy input for another 
integrator, the resulting digital curve will tend to 
straddle the average continuous value of z. Atendency 
to accumulate round-off error is reduced. 
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PART III - WHAT INTEGRATORS CAN DO 



Problems solvable with the DA- 1 will take the form of 
equations. The equations may be algebraic, differen- 
tial or integral inform, ormayoccur as simultaneous 
sets of such equations. The first step in obtaining the 
solution is to arrange the integrators to generate all 
of the terms appearing in the equations. The genera- 
tion of terms by integrators will be considered in this 
section. 

The basic integrator equations (Equations 1 and 2) can 
be represented schematically. Eachpointed block rep- 
resents one integrator. Each rectangular block repre- 
sents the constant multiplier associated with the inte- 
grator. 
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Figure 6 

Assume that prior to the receipt of any dy or dx inputs, 
the Y registers of both integrators were empty. In 
the first integrator, dy inputs accumulate to form y, 
A dz output is generated equal to ydx. It is modified 
by the constant multiplier associated with the integra- 
tor to form Kl dz. The Kl dz output is sent as the 
secondary input to another integrator where it is ac- 
cumulated in the Y register of the second integrator. 
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Now if at the start of computation the Y register of 
the first integrator had held the value Cj, and if the 
Y register of the second integrator had held a value 
C 2 , the output from the first integrator would be 

Kdz = Ki(y+Ci) dx 
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c e* K i/Y dx * K » c » x 
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Figure 7 

and there would be accumulated in the second inte- 
grator: 

C 2 + Kiz = Ki/ydx + KiC^x + C 2 . 

If the two inputs to an integrator are tied together as 
in Figure 8: 

Kdz = K(x+Ci) dx 



x+c. 




Kix+c^dx+Kdz. 



Figure 8 

which when accumulated will produce the quadratic 
function 

C 2 + Kz = K(x 2 /2+Cix) + C 2 

In the equation and diagrams in the remainder of this 
section, the constants will be omitted. These con- 
stants, constants of integration, are determined by 
the initial settings of the integrators. 
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If the output of an integrator is fed back as its own dy 
input, as in Figure 9, an exponential function is gene- 
rated: 

dz = zdx 



**nd<Z* 



Figure 9 

More complicated functions can be formed by causing 
the Y register of an integrator to accept more than one 
input. An example is shown in Figure 10. 
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Figure 10 



Mapping 



The diagrams in both Figures 6 and 7 show pairs of 
integrators with connections representing information 
flow drawn between them. The preparation of such a 
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flow diagram to represent the generation of afunction 
or the solution of a problem is called "Mapping" . 
Some examples of maps of functions follow. 

The map of a complex problem may require a great 
many integrators. It is in mapping that the creative 
aspects of programming for the DA- 1 lie. The solu- 
tion of a problemmay usuallybe represented by more 
than one map, each map representing an alternate 
approach. The procedure for mapping a problem is 
described in detail in a subsequent section. 



Generation The functions mappedon the following pages are com- 

of Functions mon ones useful in the solution of problems. The 

types which may be generated internally with Acces- 
sory DA-1 include algebraic, trigonometric, exponen- 
tial, logarithmic, hyperbolic, Bess el and probability 
functions. 

Maps for function generation can be devised by use of 
a standard table of integrals. For example, the flow 
diagram in Figure 11 is based on the relationships: 

In v - / — dv 

I = - f- d In v 
v J v 

The first equation is representedby the first integra- 
tor in the map: 

dy = d(±) 

dx = dv 

Kdz = d In v 

K = +1 

The second equation is represented by the second in- 
tegrator: 

dy = d(i) 

dx = d In v 

Kdz = d{-) 
K = -l 
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In Figure 12, the term — is generated by the same 
means as in Figure 11. In the third integrator of the 
example, the relationship is used: 



n 



j n n j i 

dv = nv d In v 



dz = dv 
dy ■ = dv 11 
dx = nd In v 



We can show this equation to be a form of de =e du: 





n n In v 

v = e 


let 


n In v = u 


then 


V n _ e u 


and 


In v = — 
n 




d In v = d - 
n 



By definition of a 
natural logarithm 



n 



And by substitution in the equation for dv : 

jz u k ne du u , 
d(e ) = = e du 

v ' n 

An interesting exercise for the reader would be to 
demonstrate the validity of the flow diagram in Fig- 
ure 13 by similar means. 
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Modes 



An integrator may be coded to have either of two 
modes of operation. The "normal" mode istheone 
which has already been described: from the dy anddx 
inputs an output dz is generated so that 

dz = ydx. 

dz is equal to either +1 increment, 0, or -lincremeit. 
The output occurs in an iteration in which the sum of 
ydx products exceeds one increment of scaled output 
value. Since y may be positive or negative, and dx 
may be positive or negative, the sum of ydx products 
may be either positive or negative. The R register 
holds the sum of the ydx products and the overflow 
of the R register provides dz. 



Decision 
Operation 



In the second mode of operation, the dz output of an 
integrator is determined differently. The R register 
is no longer used. Whenever the Y register holds a 
non-zero value and a dx input occurs, a dz output 
will occur. The sign of this output will be algebrai- 
cally determined by the product of the signs of y and 
dx. Symbolically: 

If y £ 0: dz = (sign of y) (dx) 
If y = 0: dz = 0. 

This mode is called "Decision" operation. It is use- 
ful in simulating the action of a servo, for the addition 
of functions, and in handling discontinuities in func- 
tions. 



Servo 
Operation 



One use of a decision integrator is called "Servo" 
operation. An integrator coded as a decision element 
can be made to act as a servo if it is also codedto 
receive a dx input every iteration. A servo wouldbe 
used to continually make a correction that tends to 
reduce an error, 

A servo may be used to generate u as a function of v 
where u is defined implicitly as a function of v by 
F(u, v) =0. For example in Figure 19, a square rod! 
is generated by the determination of u in the relation- 
ship 

v - u 2 = 0. 
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In differential form, this equation becomes 

dv - 2u du = 0. 

dv and -2u du are made the dy input to the integrator 
coded as a servo. 

When incremental changes in v cause F(u, v) ^ 0, the 
servo generates incremental changes in u until 
F(u,v) = 0. B 
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do 






u =.-Vv 



Figure 19 

Figure 20 shows a hook-up for the case in which 
F(u, v) = e u -v. 
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Figure 20 
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Figure 21 shows the case when 

F(u, v) = y(v) cos u - x(v) sin u. 

-1 y 
In this latter case u = tan (£) and 

^< u ' v > <0if Ixl +|y|>0. 
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Figure 21 
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Addition of If dz's from different sources are to be combined as 

Functions the dy input to a single integrator, the various dz's 

may be connected directly to the dy input line. The 
only stipulation is that the total value of the dy input 
during a single iteration must be between the range of 
-7 to +7 increments. 

If dz's from different sources are to be combined as 
the dx input to an integrator, they first must be passed 
through a decision integrator in the manner illustrated 
in Figure 22. When connected in this manner, the in- 
tegrator is called an "adder". The Kdz output of the 
adder is the negative of the sum of the inputs. 
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Figure 22 
"Adder" Connection for Addition of Functions 

Figure 23 shows an adder used to obtain du = +dx+dw. 
As a result of the servo-like operation of an adder, 
the function u will be such that x + w - u = 0. 
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Figure 23 
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Disconti- A decision integrator may also be used to generate a 

nuities wide variety of discontinuous and non-linear functions, 

such as square waves, saw-tooth waves, etc. 

Figures 24 and 25 show two integrator hook-ups that 
involve decision integrators. In Figure 24 a decision 
integrator is being employed to generate the absolute 
value of a function, u. Since, in any cycle, the same 
increment is used both as a primary and a secondary 
input, the decision integrator emits the change in the 
absolute value of u as its incremental output. In Fig- 
ure 25 two decision integrators are being used to gen- 
erate the saw-tooth functions, w and v, by means of 
same properties. 
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Figure 24 
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Figure 25 



Empirical Empirical functions may be entered into computations 

Functions in one of two ways. An empirical curve may some- 

times be approximated by the combination of known 
functions which may be generated on the DA-1. Or, 
if the empirical data is in numerical form, the inte- 
grand values may be sent to the DA-1 via the General 
Purpose Computer as it is required. 



Initial 
Conditions 



A feature of the DA-1 is its ability to remember the 
initial conditions of the problem being solved. The 
initial values are determined by the programmer from 
that particular solution of the equations which corres- 
ponds to a desired starting point for computation. (For 
example, in Figure 3 the point y , x Q specifies the ini- 
tial conditions.) The initial values, which are stored 
in the Y registeisof the appropriate integrators at the 
beginning of computation, are destroyed by the compu- 
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tation as the Y register is incremented and decre- 
mented. The initial values may be stored in a second 
set of registers called the Initial Conditions Registers . 
The contents of these registers do not change during 
the solution of the problem. At any time in the tabu- 
lation of the solution, the computer may be halted and 
initial conditions reinserted into all integrators. The 
problem may then be recomputed, with the initial con- 
ditions modified, if desired. 
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PART IV - HOW TO MAP SOLUTIONS 



Mapping 



There are three steps involved in programming a dif- 
ferential equation for the DA- 1 . The three steps are 
mapping, scaling and coding. "Mapping" consists of 
arranging an integrator information flow diagram, 
such as those in the preceding section, in a form that 
corresponds to the problem to be solved. 




<>Kd* 
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Figure 26 



dz may be sent only to the 
four consecutive following 
integrators. 



Kdz may be sent to all in- 
tegrators. 



As an illustration of mapping, consider the example 
discussed in Part I. 

After the equation is re-arranged and put in differen- 
tial form (page 8), we have 

d(^l) = -. 297 d Aiv 2 dv + d /u v du + d cos u (Eq. 3) 

dv 
If ctu is assumed to exist within an integrator, the de- 
pendent variable dv maybe obtained from the output of 
the integrator: 







+ I 



■>■ dv 



dv = j— du 
du 
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By sending dv to a second integrator as both the dx and 
dy input, d(v 2 ) may be generated: 



M * 



dv 



+ 2 



+ <d(v 2 ) 



d(v ) = 2v dv 

By similar means, d(u ) may be obtained from du: 



u 




<t(o*) 



d(u 2 ) = 2u du 

The input terms are now available to generate d^v dv 
and d/u 2 du which in turn are necessary to generate 
the first two right-hand terms of Equation 3. 

dV 




^A/a^u 



'/ 



d / u 2 du = u 2 du 



Now if each of these latter terms is used as the dx in- 
put to an integrator, and if the du and dv are used as 



respective dy inputs: , 



O 




-X9 7 



.Z97 



d/u^d' 
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7 d /uv 2 dv = -.297 uv 2 du = -.297 ud/v 2 dv 
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d/u^Vdu 



u 



vdu = u^v du = vd 



u du 



The third right-hand term of Equation 3, cos u, is an 
example of a function whichmay be generated directly 
in a differential analyzer. Two integrators are re- 
quired for the generation by making a closed loop of 
the relationships. 

d (cos u) = -sin u du 
d (sin u) = cos u du 

Now d(^J) may be formed by combination of the other 
three terms in Equation 3. d(§g) is made the depei^ 
dent differential input into the integrator in which ^. 
was assumed to exist in order to form dv. 

The completed flow diagram is shown in Figure 27. 



As a second example, consider the second order non- 
linear differential equation: 

dx 



Three flow diagrams, each representing a different 
method of solution, are shown for this problem in Fig- 
ures 28, 29 and 30. 
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dv 
du 



U 



U 



SIN U 



COS U 




Figure 27 
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Equation 



Explanation: 



25 



28 



31 



34 
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dx 2 



-P- + 

dx 



METHOD I 
2 



+ sin y + A 



The highest derivative has been separated from the 

d^y 
rest of the equation. Assume — ~ in Integrator 25 in 

dv ^ dv 

order to form d(-pp-). Accumulate -^- in Integrator 28 

and generate dy. Accumulate y in Integrator 31 and 
form d(y^); Integrators 34 and 35 generate d(sin y). 
The terms are added together and fed back into Inte- 
grator 25 to close the loop. 
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Figure 21 
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METHOD II 



Equation 
Explanation: 



d y + dy _ 2 _ 



dx< 



dx 



sin y — A = 



A decision integrator is used as a servo to control the 
highest derivative in such a manner as to hold the sum 
of the terms within one digit of zero. Comparative 
accuracy is slightly less than that obtained by the first 
method. This method has the advantage that the equa- 
tion need not be explicitly solved for the highest deri- 
vative. K , 




Figure 29 
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METHOD III 



Equation: 



i?" 



$L + y 2 + sin y + A 
dx 



Explanation: 



If we integrate with respect to x and then obtain the 
differential (as in the example in Part I), the terms 
may all be generated by integrators: 



*&• 



dy + (y + sin y + A) dx 



d^y . 
This method has the possible disadvantage that j is 

dx 
not explicitly generated in the computation. However, 

the method layout often conserves integrators. 
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Figure 30 
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A similar procedure may be used in mapping a prob- 
lem involving a simultaneous set of differential equa- 
tions, as shown in Figure 31. 



d w 



dt" 



dv 



3 
dv =w^ + vt 
d?^ dt 




d (d_£ ) = wdw + vtd+ 

at ' 



dv + vdt 



Figure 31 
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It may happen that a particular scheme for mapping 
will give rise to an integrand which becomes infinite 
at some point in the problem. For the equation, 



d 2 y 



_ l_ dy _ 



dx 



the map shown in Figure 32 cannot be used when x is 
small or zero. Figure 33 shows an alternative map 
that may be used in the region when x is small, du is 
defined as * dy in the problem and is used in a manner 
which prevents an integrand from increasing beyond 
bounds. 
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Figure 32 
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Figure 33 

A simple map may often be found for a problem by per- 
forming a substitution of variables. For example, in 

the equation, 






ln x, 



the machine independent variable may be us ed for d In x 
as shown in Figure 34. 
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dt"d In x 




Figure 34 

Note that machine time, dt, is not always used as an 
independent variable in the problem. Machine time 
may even be used as the problem dependent variable 
as shown in Figure 35. 

The map shown in Figure 35 can be used with the ini- 
tial condition, x = 1, y = 0, even though 

^dx* 

is infinite. 
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Multiplication 
of Variables 



In most cases the product required can be placed in 
differential form; products may then be built up using 
just one integrator per multiplication. For example, 
the term uvywdx may be generated with four integra- 
tors as shown in Figure 36. 



W 




^nv^wdx 



Figure 36 



In cases in which the product cannot be placed in dif- 
ferential form, use is made of the relation 



uv 



udv 



= /vdu + j\ 
and the corresponding map is shown in Figure 37. 



42 



ov 




Figure 37 



In Figure 1 at the front of this manual, a plot is shown 

clx 
of an elliptical spiral. The picture shows ^- plotted 

against x in the equation 

l^x + dx +o 2 x = 

Ht 1 w 

The subsequent map is shown in Figure 38. 
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Figure 38 
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PART V - HOW TO SCALE NUMERICAL VALUES 

There are two additional steps to the preparation of a 
program for the DA-1 after the flow diagram has been 
prepared. They are the selection of scaling exponents 
for each integrator and the coding (putting in computer 
language) of the information from the scaled flow dia- 
gram. Scaling will be discussed in this section. 

Choice of Three binary scale exponents must be specified for 

Srale every integrator. One, S v , is determinedby the man- 

powers imum value expected to be reached by the integrand 

v; one, Sv, is determined by the value of the constant 
multiplier, K; the third, S dy , specifies the magnitude 
of variation iny to be represented by a dy increment. 

Consider in turn the method of obtaining each: 

Sv TofindS v , the largest magnitude expectedto be reached 

7 by the integrand is estimated. S y is the power of 2 

such that 2^ is larger than this maximum. (If "* com- 
putation the integrand tries to reach or exceed 2 y, 
the DA-1 will automatic ally halt and the overflow neon 
on the main computer will light. ) 

Example: The y variable is the speed of an auto- 
mobile which is known to remain less than 100 
miles per hour. 

Since 2 6 = 64 and 2 7 = 128, the smallest scale 
exponent that could be used for S y is 7. 

S k is the power of 2 such that 2 S k is larger than the 
constant multiplier, K. 

Ex ample: If K = 296, find S k . Since 2 8 = 256 
and 2* = 512, S k is set equal to 9. 

Sh v represents the increment of the Y variable. It 
sSLifies the value of a single dy input, the smallest 
variation to be permitted in the integrand. S dy is 
chosen so that 2 s dy is the nearest power of 2 to the 
desired smallest variation. 

Example: In the automobile problem above, the 
speed is required to the nearest tenth of a mile 
or better. One -tenth is between one -eighth and 
one -sixteenth. Since 1/16 mph = 2"4 mph, set 
Sdy to be equal to -4. 

If an integrator is used as a constant multiplier, make 
S dy = S y - 26. 



Sk 



S dy 
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All exponents in DA-1 programming are expressed in 
excess -fifty form. In the preceding three examples, 
therefore, Sy would be written as 57, Sfc as 59 and S<j v 



as 



46, 



Relationship 
Between 
Scaling 
Exponents 
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Figure 39 



The first step in scaling a problem is to estimate Sy, 
S^y and S^. for each integrator -multiplier pair* It 
is then necessary to determine the compatibility of 
the scale exponents estimated for each integrator with 
those of the other integrators. To do so, the addi- 
tional scaling exponents listed in Figure 39 must be 
calculated. The scaling exponent which exists at the 
output of an integrator (S z \ or S z i) must be equal to 
the S^y or S x for the integrator to which it goes as an 
input, S z i and S z 2 are determined by the formulas in 
Figure 39. 

It may prove necessary to modify some of the estimated 
scaling exponents in order to obtain compatibility. 

An example of the procedure for the adjustment of scal- 
ing exponents is shown in Figure 40, In Integrator B, 
S^y was first estimatedto be -4; however, the calcula- 
tion of S z 2 for Integrator A necessitated that S^y for 
Integrator B be changed to -6. 
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Figure 40 

Figures 41 and 42 show the same flow diagram as Fig- 
ure 40 with the necessary integrator scaling adjust- 
ment made indifferent manners. In Figure 41 the scal- 
ing exponent of Y a in Integrator A was the one changed 
to obtain compatibility. In Figure 42 the scaling ex- 
ponent of K a of Integrator A was changedto accomplish 
the same result. 



IVa^2 u 



Q> 



Ye> 




2^=2 



-4 



Figure 41 
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Figure 42 

When several integrators provide dy inputs to another 
integrator, the output scaling exponents of each of the 
sources must be equal to Sdyfor the destination inte- 
grator. 

The scaling exponent of the machine independent varia- 
ble, dt, may be arbitrarily chosen in order to make 
the other scaling exponents realizable and practical. 
A convenient preliminary value for the scaling expo 
nent of dt is the average of the values of S^y 
the integrators in the program. 



S y for 



The scaling relationships described are summarized 
in Figure 43. 
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Note: These scaling relationships do not apply to inte- 
grators coded for decision element operation. 



Figure 43 
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PART VI - HOW TO CODE AND OBTAIN 
NUMERICAL SOLUTIONS 

After the problem has beenmappedby means of a flow 
diagram, after an initial value for y and a value for k 
has been chosen for each integrator, and after scaling 
has been completed, the problem may be coded. In 
coding the scaled flow diagram is transcribed into a 
set of numerical codes understandable by the computer. 

Rules and considerations for the assignment of inte- 
grator numbers and for other aspects of coding are 
given in Table I. Each integrator in the problem is 
coded by use of the form in Table II. Read-out instruc- 
tions to the computer are coded in the form given in 
Table III. The information from the coding forms may 
then be entered into the computer by following the pro- 
cedure in Table IV. 

If desired, the program may be modified after being 
filled into the computer. Operation codes with which 
to do so and other codes for instruction of the DA-1 
are listed in Table V. 

Since it is quite easy to violate one of the coding rules 
and procedures while familiarity is being gained with 
them, facilities are included for the automatic detec- 
tion of errors. 

An error during typein of the program will cause an 
immediate typeout of a number indicative of the type 
of error. In Table VI the possible errors and their 
indicator numbers are tabulated. 
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TABLE I - Coding Considerations and Rules 



1. Integrator 00 automatically has a dx input every 
iteration. The output of Integrator 00 (dt) is the 
independent variable that has been called "machine 
time". 

2 6 If Integrator 00 receives a programmed dx input, 
positive or negative, computation in theDA-1 will 
automatically halt. Therefore, the output of an 
integrator maybe used to halt computation by cod- 
ing the integrator to be the dx input to Integrator 
00. 

3. An integrator may be coded to receive a dx input 
from one other integrator. 

4. An integrator may be coded to receive dy inputs 
from any number of integrators. The total value 
of dy inputs in any one iteration must lie between 
-7 and +7 increments. Therefore, if information 
is not known about the total value of dy inputs to 
an integrator during an iteration, no more than 
seven integrators should be coded for dy input. 

5. The kdz output of an integrator (see Figure 44 ) 
may be sent to the input line of any integrator and 
any number of integrators. The kdz output sent 
to the integrators, numbered one greater and two 
greater than the integrator providing the output, 
will be one iteration old. 

6. The dz output of an integrator (see Figure 44 ) 
maybe sent to the integrators numberedone great- 
er, two greater, three greater, and four greater 
than the integrator providing the output. 



M 




Maybe sent to any 
^ YAi integrator, includ- 
ing N. Output from 
N to Integrators 
N+l and N+2 will 
be one iteration old. 



May be sent to 
Integrators N+l, 
N+2, N+3 and 
N+4. 

Figure 44 
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TABLE I - Coding Considerations and Rules (cont, ) 



7. Numerical values for y and K are expressed in 
floating decimal point form. The numerical por- 
tion is written as decimal point followed by seven 
digits. The exponent portion may range between 
10" 4 and 10+9. The exponent of 10 is expressed 
as an excess -fifty number preceding the numerical 
portion. 

Example: The number "4963 " 
~~ wouldbe expressed in 

computer language as 
"54 4963000" 

8. The DA-1 maybe instructed to operate as a 27 in- 
tegrator machine, a 54 integrator machine, an 81 
integrator machine, or a 108 integrator machine. 
The integrator numbers which may be used in each 
case are: 

27 integrator operation: Integrators 0, 4, 8, 

12 ... . 104 

54 integrator operation: Integrators listed 

above, plus 1, 5, 9, 

13 ... . 105 

81 integrator operation: Integrators listed 

above, plus 2, 6, 

10, 14 ... . 106 
108 integrator operation: Integrators listed 

above, plus 3, 7, 

11, 15 ... . 107 

9. The numerical value of the difference between 
the scaling exponents SyandSdy (that is, Sy - S<jy) 
must be between 3 and 26 inclusive. 

10. The decimal exponent of an integrand coded for 
typeout must be large enough to permit typing the 
largest possible integrand value as a fraction. 

Example: If the binary scale ex- 
ponent, Sy, equals 10, 
the integrand could be- 
come larger than 1000 
but will remain under 
10,000. Therefore the 
smallest permissible 
value of the decimal ex- 
ponent of the integrand 
would be +4. Larger 
values could be used. 
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TABLE II - Coding Form 

± 50 + M 



_L_ 



± 50 + NI 



5Q + 5y 50 + S 
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Yq and K are expressed as sign, decimal exponent in 
excess-fifty form, and decimal fraction. The number 
"+8263", for example, would be written" +54.8263000". 

"I" is the number of the integrator being filled. 

"dx" is the number of the integrator from which dx 

comes. 

"dy" is the number of an integrator from which dy 
comes • Though space for only seven dy inputs are 
shown, more may be filled if desired. 

The "Decision" line should not be typed unless the in- 
tegrator is to be a decision integrator. 

Integrator numbers are always written with two digits. 
Integrator numbers 100 to 107 are written as uO to u7. 



*The earlier model typewriter has no (s) key; 
wherever the key is indicated, use the 
"s" key on the alphabetic keyboard. 



52 



TABLE II - Coding Form (cont. ) 

Adxordy input specified as 00 to 99 or as uO to u7, 
refers to the multiplied (kdz) output of the integrator 
from which the input comes. To specify that the input 
to any integrator, N, be from anunmultiplied (dz) in- 
tegrator output, write the integrator number of the dx 
or dy input as follows: 

Write zl if input from Integrator N-l 
Write z2 if input from Integrator N-2 
Write z3 if input from Integrator N-3 
Write z4 if input from Integrator N-4. 
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TABLE III - Readout Control 



1. To list an integrator for typeout 
SO+*i 50-*- Sy NIJ MgEg 



g 



TAB 



The contents of the Y registers of any twelve, 
or fewer, integrators may be automatically 
typed out at selected intervals, 

Typeout is decimal, and is expressedin float- 
ing point form as a decimal fraction. 



2. To code typeout interval 



TVREOUT 

CONTROL. 
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Any integrator may be selected to be typeout 
control integrator. Integrands listedfor type- 
out will be typed out whenever the typeout 
control integrator has a dz output. 

Explanation: As stated in Table IV,computa- 
tion begins with the command "x-^TAB)®''. 
This command causes typeout followed by re- 
start of computation when the DA-1 is halted. 
In the first of the two coding steps above, the 
R register of the typeout control integrator is 
set to zero (rather than . 5), so that all type- 
out intervals will be equal; in the second, 
Integrator 00 is coded to have the typeout 
control integrator as its dx input. This second 
step causes computation to halt when the type- 
out control integrator has an output. 
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TABLE IV - Operating Instructions 

1. Initially, the three switches on the typewriter base 
("Compute", "Enable" and "Punch") should be off. 

2. Insert magazine DAFPER-1A in rewound position 
on the photo reader, 

3. Plug in Accessory DA-1 (two plugs). Turn on 
power switch on back panel of DA-1. 

4. Set switches on back panel of DA-1 as follows: 

Switches 9 and 13 on, others off, for 27 inte- 
grator operation* 

Switches 8, 9, 12 and 13 on, others off, for 54 
integrator operation. 

Switches 7, 8, 9, 11, 12 and 13 on, others off, 
for 81 integrator operation. 

All switches on for 108 integrator operation. 

5. TurnonG-15 Start Switch. Wait several minutes 
for warm-up cycle to be completed. 

6. Hold RESET button in until the DC voltage is turned 
on. 

Two blocks of tape will be read in. The first 
block is the Number Track for the G-15 Com- 
puter. The second block is a loading and mem- 
ory - setup routine. Initial values of zero are 
put in Y registers and values of +1 are put in 
K registers. 

7. Wait until the green "Ready" light turns on. Then 
put Compute switch on typewriter base to GO. 

The typewriter will execute a carriage return. 
Neons on G-15 power control panel will be in 
TYPEIN configuration. 

8. Type "N TAB(§y ! where N refers to the number of 
integrators in one iteration according to the fol- 
lowing code: 

N=l for 27 integrator operation 
N=2 for 54 integrator operation 
N=3 for 81 integrator operation 
N=4 for 108 integrator operation 

The computer will type outthe number of inte- 
grators and additional tape will be read into 
the computer . When tape reading has been com- 
pleted, the typewriter will execute a carriage 
return. Neons on the G-15 power control panel 
will be in the TYPEIN configuration. 
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TABLE IV - Operating Instructions (cont. ) 

9. The coded program may now be typed directly from 
the coding forms of Tables II and III. All contents 
of the forms which are in blocks must be typed. 
The remaining information is merely explanatory. 

Input information is automatically converted 
into proper binary form and storedon the mem- 
ory drum. Initial values of integrands are 
stored in Initial Conditions Registers on the 
drum. After the filling of integrand or constant 
multiplier values, the computer will retype the 
number in hexadecimal notation. After any fill 
operation, the computer will retype the opera- 
tion code. 

10. Type "y-KTAB)®" 

This instruction inserts initial values of all in- 
tegrands from the Initial Conditions Registers 
into the Y registers. 

1 1 . For reference during readout, the integrator num- 
bers and exponents of integrators to be read out 
may be automatically typed by the computer. 

Type "u-(TAB)<§r 

Integrator numbers of integrators to be 
read out will be typed out. 

Type M v-{TAB)(3) M 

Binary scale exponents of read-out in- 
tegrators will be typed out. 

Type "v+tTAB)® 1 

Decimal exponents of read-out integra- 
tors will be typed out. 

12. Type "x+tTAB)® 1 

The computer will type out initial values of in- 
tegrands and will then start computing, typing 
out computed results in decimal form at the 
programmed intervals. 

Halt Control 

Computation will halt automatically either when Inte- 
grator 00 receives a programmed dx input or an inte- 
grand overflow occurs. If it is desired to halt compu- 
tation at other times, and thus permit typein of further 
instructions: 
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TABLE IV - Operating Instructions (contc ) 

1 Move Compute switch to BP when control -panel 
neons are in the TEST DA-1 OFF or TEST 
READY configurations. 

The DA-1 will halt after the next typeout. 

2. After the typeout is completed, move the Com- 
pute switch to OFF. 

3. With the Enable switch on, type "f". 

4. Move the Compute switch to GO, 

After one more typeout, the control- 
panel neons will revert to the TYPEIN 
configuration. 

To halt computation in case of program difficulties: 

1. With the Compute switch OFF and the Enable 
switch on, type w @cf ". 

2. Move the Compute switch to GO, 

The control- panel neons will be in the 
TYPEIN configuration. 



Memory Punchout 

The contents of DA-1 memory can be punched on tape. 
A problem in memory may be punched out at any stage 
of computation and later reloaded. Computation may 
then proceed either from the beginning of the prob- 
lem or from the point where it was interrupted* If 
computation is to be resumed exactly from the point 
at which it was interrupted, dx and dy inputs should 
not be addressed to integrators 00, 01, 02 and 03. 

The procedure for initiating memory punchout is: 

When the neons on the power control panel are in 
the TYPEIN configuration (follow Steps 1-4 under 
Halt Control , above, to obtain this configuration), 
type "9 +1TAB)®". 
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TABLE IV - Operating Instructions (cont. ) 



The contents of DA-1 memory will be 
punched on tape. A bell will ring after 
about ten minutes to indicate the last 
block of tape being punched. 

Notes: Initial conditions are automatically re- 

inserted in the integrators after punching 
is completed. Therefore, the problem 
is in its original state in the memory of 
the computer. 



Reloading Punched Out Contents of Memory 

The contents of DA-1 memory that were punched on 
tape can be re-entered into the computer. After 
DAPPER-lAhas been read in, the procedure is: 

1. Move the Compute switch to OFF and then re- 
move Magazine D APP ER- 1 A without r ewinding . 

2. Insert the tape containing the punched out con- 
tents of memory on the photo-reader. 

3. With the Enable switch on, type "p". 

Wait for the photo-reader light to gooff. 

If the bell rings and the computer halts, 
a reading error has occurred. In this 
case, rewind the tape and repeat Step 3, 

4. Move the Compute switch to GO. 

The subsequent blocks on the tape are 
read in. In case of error the tape is 
automatically reversed, the bell rings, 
and the block is re-read. 
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TABLE IV - Operating Instructions (cont. ) 

5 U Check the switch settings on the back panel of 
the DA-1 to be sure the switches are in the 
same positions they were before the memory 
punchout of the problem. (See TABLE IV - 
Operating Instructions, Step 4, page 55. ) 

6. Replace Magazine DAPPER-1 A on the photo- 
reader . 

Input information may be typed under 
control of the DAPPER-1 A routine. 

7. To insert initial conditions, type "y(TAB)® 11 . 

8. To resume DA-1 computation, from the point 
where memory was punched out, omit Step 7 

and type "x(TAB)@". 



Monitor Operation 

Programs may be monitored if desired. In this oper- 
ation the contents of readout integrator s are automati- 
cally typed out every iteration. The procedure is: 

1. Move Compute switch to GO. 

2. Type "u +tTAB)<§>". 

The current values of readout integrands 
will be typed. On completion of the type - 
out, the control -panel neons will revert 
to the TYPEIN configuration, 

3. Type "x -(TAB)®". 

The program will now be executed and 
readout integrands will be typed out each 
iteration. 

To stop monitoring, the procedure is: 

1. Move the Compute switch to BP when the con- 
trol-panel neons are in the TEST READY con- 
figuration. 

The computer will execute one more 
type out and then halt. 
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TABLE IV - Operating Instructions (cont. ) 

2. When the computer has halted, move Compute 
switch to OFF o 

3. With the Enable switch on, type M <§)cf". 

The neons will revert to the T YPEIN con- 
figuration. 

To skip typeout intervals while monitoring: 

In order to monitor program execution at points 
away from the start of computation, move Com- 
pute switch to BP and type "x-lTAB)®! 1 . Moving 
the Compute switch to OFF and back to BP will 
alternately start and stop the DA-1. At each halt, 
readout integrands will be typed out. Normal mon- 
itor operation can be resumedby moving the Com- 
pute switch to GO o 

Caution: Move the Compute switch only when the 
Halt neon is lit or control-panel neons 
are in the TEST READY configuration. 



Overflow Detection 

After M x-KTAB)(D M has been typed: 

If an integrand overflows during normal computa- 
tion, the DA-1 will halt. The photo-reader will 
read in a section of tape that contains overflow in- 
structions. The code digits "zz" will be typed out 
to indicate that the halt was due to an overflow. 
The integrand which overflowed will be typed out 
in hexadecimal notation. The integrator number 
in which the overflow occurred will be typed out in 
decimal notation. The integrand which overflowed 
will be cleared to zero. The current values of all 
readout integrands will be typed out. 

After n x-(TAB)® M has been typed: 

If an integrand overflows during monitor operation, 
the DA-1 will halt and the Overflow neon will light. 
However, in order for the overflow detection and 
typeout operations listed above to occur, "T^TAE)®" 
must now be typed. 



60 



I 






3 - tab® 


1 




4 -! T4B © 


1 


— — 

1 

i 


8|-H® 



TABLE V - Operation Codes for Program Control 
Remove dx input from Integrator I. 
Remove dy inputs from Integrator I. 



Change Integrator I from decision use to 
normal use. 



JAB [5 



7|+ T*8|© 



1?1 + M£ 



u 



^s 



Remove list of integrators coded for read- 
out. 

Type out integrand (hexadecimal nota- 
tion) and integrator number (decimal no- 
tation) which caused overflow. Clear 
overflow integrand to zero. Type out in- 
tegrands of readout integrators, Turn off 
overflow neon. 

Punch out contents of DA-1 memory on tape. 

Type out decimal portions cf readout inte- 
grands. 



U 



TAB 



m 



Type out integrator numbers of readout 
integrators. 



V 



+ TA6 



m 



Type out tens exponent of readout inte- 
grands in excess -fifty form. 



- tab rs 



Type out Sy of readout integrands in ex- 
cess-fifty form. 



W 



TAB (5) 



Recall DAPPER-1A loading routine, 



X + 746 (S 



Type out decimal portions c£ readout inte- 
grands and compute, typing out at pro- 
grammed intervals. 



X - TAB 



Compute, typing out decimal portions of 
readout integrands every iteration. 



H© 



Insert initial values of integrands in all 
integrators. 
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TABLE VI - Error Indicators During Typein 



Typed-Out 
Indicator 



59 



58 
or 
57 



56 



53 



54 
or 
53 



52 



51 



49 



zzzzzzz 



zOOOOOO 

or 
GOOOOOz 



No typeout 

after 

M u(TAB)<D n 
or 

"xlTABldy 1 



Operation 
in Error 



1+ Fill Y 
2+ Fill K 



1+ Fill Y 
2+ Fill K 



1+ Fill Y f 



dx or dy (not 
coded 1+, 2+, 
3+ or 4+) 



3 + 

4+ 



3+ 
4+ 



1 + 
2+ 



3 + 
4+ 



u+ 
x+ 
x- 



u+ 
x+ 
z- 



dx input 
dy input 



dx input 
dy input 



Fill Y 
Fill K 



dx input 
dy input 



u+ Typeout 
x+ Typeout 
z+ Typeout 



Type of Error 



(50+N) exceeds limits 



Sy too large or too small 
or scaling in error 



S y - S^y beyond limits 



Integrator number 
greater than 107 



Integrator number of 
source not legitimate 
numbe r . 



Integrator number of 
destination greater 
than 107 



Integrator number 
greater than 108 



Machine error 



Overflow in computation 
preceding readout 



Sy in readout list is too 
large or too small 



(50+N) in readout list 
exceeds limits 
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TABLE VI - Error Indicators During Typein (cont. ) 



The decimal exponent of a readout integrand must be 
large enough to permit typing the largest possible in- 
tegrand value as a fraction, (See Table I.) If 50+N is 
too small or much too large, M z000000 M or "OOOOOOz" 
would be the typed out error indicator. 

If an error indicator is typed out accompanied by two 
carriage returns, rather than one, the type of error 
is not determinedby the indicator number. The error 
is one in which a zero preceded '{TAB)® 11 . Zero in that 
position corresponds to no legitimate operation code. 
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TABLE VII - Graphical Output 

Facilities are included intheDA-1 for the use of graph- 
plotter output. 

A graph-plotter accessory is available which receives 
Ay and Ax increments from the DA-1 and plots them 
in 0. 01 inch increments on a standard roll of paper, 1 
foot by 100 feet. 

Programming specifications for the use of the plotter 
are summarized below. 

Output to Graph-Plotter Plugged Into Graph -Plotter 
Jack 

^y The output of any integrator in the series: 

1, 5, 9, 13, 17, .... 105 

may provide the Ay output to the plotter. The 
integrator will provide Ay if it is coded to be 
the dx input into Integrator 03. 

Z^x The output of any integrator of the two series: 

0, 4, 8, 12, 16, 104 

1, 5, 9, 13, 17, .... 105 

may provide the Ax output to the plotter. The 
integrator will provide Ax if it is coded to be 
the dx input into Integrator 02. 



Output to Graph-Plotter Plugged Into Graph -Plotter- 
Follower Jack 

Ay The output of any integrator may provide the 

£y output, except integrators in the series: 

2, 6, 10, 14, 18, 106. 

The integrator will provide Ay if it is coded to 
be the dx input into Integrator 01. 

/\x The unmultiplied output of Integrator 107 will 

provide the Ax output to the plotter. No addi- 
tional coding is necessary. 



64 



TABLE VIII - Control Panel Neon Configurations 



Significance 


Pattern 


TYPEIN 


O O Characteristic 
O O O Command Line 
# O O Source 
9 Destination 
• O Input -Output 


TEST READY 


O O Characteristic 

O O Source 

Destination 


TEST DA-1 OFF 


Characteristic 
O O Command Line 
O O Source 
Destination 


TEST PUNCH 
SWITCH 


O w Characteristic 

O O • Source 

Destination 


START DA-1 


O O Characteristic 
O O Command Line 
O O • Source 
Destination 


TYPEOUT 


O O Characteristic 
O O Command Line 
O O Source 
Destination 
O O O • Input -Output 
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TABLE VIII - Control Panel Neon Configurations (cont„ ) 



"Typein", M Test Ready" and "Test DA-1 Off" configu- 
rations occur at the times specified in Table IV, 

The "Test Punch Switch" pattern means that the Punch 
switch on the typewriter base is on. The switch shouLd 
be turned off for Routine DAPPER- 1 A. 

The "Start DA-1" or the "Typeout" pattern will occur 
when the Compute switch is in the BP position after 
"x-|TAB)<D". 

The "Typeout" configuration will occur when the Com- 
pute switch is in the BP position after "x-KTAB)^ 1 . 
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POINTS TO REMEMBER 



1. The value in the Y register of an integrator coded 
as a servo, or as an adder, must remain close to 
zero throughout computation for results to be ac- 
curate, 

2. If the value in a K register is to be +1, that value 
need not be inserted by the programmer. Routine 
DAPPER-lAputs +1 in all K registers unless some 
other value is typed in* 

3. If a mistake occurs during the type in of dx or dy 
inputs, remove the inputs in error by using the 
proper operation code from Table V before retyping 
the correct input information. 

4. In turning on the G-15 and DA- 1 combination, be 
sure that the warm-up cycle in both units is com- 
pleted before turning on DC voltage. 

5. If the computer is used when Accessory DA-1 is 
not connected, the dummy plug must be in place 
in back of the computer. Whenever the computer 
is in use and Accessory DA-1 is connected, the 
Power switch in the DA-1 must be on. 

Computer power must be off while the DA-1 is 
being attached or unattached. 

6. To halt and return to the "Permit Typein" state 
in case of program difficulties, type M ©cf ,f with 
the Enable switch on and the Compute switch off. 



7. Routine DAPPER-lAputs initial values of zero in 
the Y Registers. There can be no output from 
Integrator 00 unless it is programmed. Accord- 
ingly, it is necessary to enter a constant, nor- 
mally + 1, in the integrand of Integrator 00 so that 
an output may occur. 
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APPENDIX 



The Program Preparation Routine, DAPPER-1A 

The scaling and coding instructions of Parts V and VI 
are based on the use of the phototape magazine DAP- 
PER-1A, Differential Analyzer Program Preparation 
and Execution Routine, Number 1A. This routine has 

been designed for typewriter input and typewriter or 
graph -plotter output. A different program preparation 
routine would provide different input-output character- 
istics . 



Internal Information Storage With Accessory DA-1 

Accessory DA-1 utilizes the memory drum of the gen- 
eral purpose computer for information storage during 
computation. Each of the 108 word times on the drum 
periphery corresponds to one integrator time. It is 
necessary to know the distribution and nature of the 
information stored on the drum in o^der to code the 
DA-1 without the use of a program preparation routine, 
or in order to write a program which simultaneously 
uses the differential analyzer and general computation 
facilities of the G-15 system. 

Form of Numeric ad information for the DA-1 is arranged dif- 

Information ferently fromnumerical information for general com- 

putation. 

The basic G-15 word consists of a sign bit, in bit po- 
sition 1 and twenty-eight bits of magnitude, in bit posi- 
tions 2 through 29. The DA-1 word consists of twenty- 
seven bits of magnitude, in bit positions 1 through 27, 
and a sign bit in bit position 28; bit position 29 is blank. 
(The presence of a "one" in bit position 29 of line 17 
during computation indicates the overflow of the regis- 
ter represented by the word,) 

In the DA-1 word, a negative number is written in 
complemented form and a negative sign is represented 
by a "one 11 bit in bit position 28, the sign position. 
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Memory There are twenty 108 -word channels or M long lines' 1 on 

Drum the magnetic drum in the G-15 Compute r. When the 

Utilization DA-1 is attached, these channels are used as follows: 

Line 

Control Information 18 

Y Registers 17 

R Registers *6 

K Registers 15 

K r Registers 14 

DY Addresses(If DAPPER-lAis used) 

27 Integrator Operation 13 

54 Integrator Operation 12, 13 

81 Integrator Operation 11,12,13 

108 Integrator Operation 10, 11, 12, 13 

DX Addresses (If DAPPER-lAis used) 

27 Integrator Operation 9 

54 Integrator Operation 8, 9 

81 Integrator Operation 7, 8, 9 

108 Integrator Operation 6, 7, 8, 9 

Y Q Registers (If DAPPER-lAis used) 

27, 54 or 81 Integrator Operation 10 
108 Integrator Operation 5 

Routine DAPPER-1A requires for its program the re- 
maining long lines on the drum. Long memory lines 
available for general use by the programmer when the 
DA-1 is in operation are summarized below, 

27 Integrator Operation 5, 6, 7, 8, 1 1, 12 

54 Integrator Operation 5, 6, 7, 11 

81 Integrator Operation 5, 6 



Line 18 Line 18 contains control pulses which are automatically 

filled by Routine DAPPER-1A. The filling of these pulses 
must be programmed if DAPPER-1A is not used. Every 
word (one word time corresponds to one integrator 
time) contains either two or three control pulses. A 
"one 11 must be inserted in bit position 26 of every word 
in the line* If an integrator is to be used as a decision 
integrator, a "one" must be placed in bit position 28 
of the correspondingly-numbered word in Line 18 In 
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each wordthat corresponds to an integrator used in the 
program, a "one" must be placed in the bit position, 
relative to the binary point, the number of which equals 
Sdy for the integrator. If the integrator is one not used 
in the program, a "one" is placed in bit position 5. 

Example: The binary point of an integrand is be- 
tween bit positions 15 and 16. If S^y = 
+4, a one is placed in bit position 19of 
the correspondingly-numbered word in 
line 18, If Sjy = -5, a one would be 
placed in bit position 10 instead. 



Line 17 Line 17 holds the Y registers or integrands. Bit po- 

sition 29 of every word is blank; bit position 28 holds 
the sign; with routine DAPPER-1A, digit positions 2 
through 27 hold the magnitude in binary form. The 
location of the binary point in the integrand is deter- 
mined by the scaling exponent Sy. S y is equal to the 
number of binary positions (not including bit positions 
28 or 29) to the left of the binary point. A negative 
number is written in complemented form. 

Example: If y = +1101. 1001 and Sy = 6, in an in- 
tegrator, the correspondingly- numbered 
word on line 17 would contain: 



OOtOOllOli. 100100000000000000000 

f \ r 

^ BIT ^ -SIX bit 

PosmON BINARY Poi.T.OK 

If y = -1101.1001 and Sy = 6, in an in- 
tegrator, the correspondingly -numbered 
word on line 17 would contain: 

01110010.011100000000000000000 

!^\ t 

bit sigtn air 

POSITION POSITION 

M 1 



Routine DAPPER-1A automatic ally performs the con- 
version from the decimal input program to the form 
which has been described. 
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Line 16 



Line 16 holds the R registers, As stated in the text, 
before DA-1 computation is begun, avalue of one-half 
in DA-1 language* shouldbe inserted in the R register 
of every integrator used in the program except the 
typeout control integrator. With Routine DAPPER- 
the insertion is automatic when the operation code 
u y(TAB)(D n is typed. 



Lines 15 
and 14 



Lines 15 and 14 hold the K and K r registers, respec- 
tively, in forms which are exactly like the forms of Y 
and R in Lines 17 and 16. A value of one -half should 
be inserted in the K r registers of integrators used in 
the program. 



The DA-1 uses "short lines" 21 and ZZ on the memory 
drum for integrator output storage , In combined oper- 
ation of the DA-1 and G-15 lines 21 and ZZ are not 
available when the DA-1 is on; if the contents of these 
lines are disturbed when the DA-1 is off, computation 
can not be resumed at the same point. 



Dy and Dx 
Addresses 



The location and insertion of pulses on the drum to 
signify the dy and dx addressing of each integrator is 
best done automatically by DAPPER-lAor some simi- 
lar routine. If it is desired to insert the pulses di- 
rectly, the required pulse positions on the drum may 
be located easily with the aid of the DA-1 addressing 
wheel. 

Pulses must be placed in bit position 29 of words N - 1 
(mod 4) in the drum channel specified by the wheel. 
The wheel, a form of circular slide rule, is available 
on request. 



Memory 
Punchout 



Lines 05 through 18 are punched on tape during a mem- 
ory punchout if a problem is set for 108 -integrator 
operation. If a problem is set for less than 108 inte- 
grators, Lines 07 through 18 are punched. 



*2000000 in hexadecimal form. 
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NUMERICAL VALUES FOR PLOT 
OF ELLIPTICAL SPIRAL 

The plot is shown in Figure 1„ 

The flow diagram is shown in Figure 38, 

The numerical values used to obtain this graphical 
solution are: 

C = .0625 

03= 1 

d 2 x 
Initial value of — -r- - . 5 

dt 

dx 
Initial value of -*— = 
dt 



S dy = -6 

AY to plotter = .5625 d(— j 

AX to plotter = . 9375 dx 
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